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The thermally induced denaturation of DNA in the presence of attractive solid surface is studied. 
The two strands of DNA are modeled via two coupled flexible chains without volume interactions. If 
the two strands are adsorbed on the surface, the denaturation phase-transition disappears. Instead, 
there is a smooth crossover to a weakly naturated state. Our second conclusion is that even when 
, the inter-strand attraction alone is too weak for creating a naturated state at the given temperature, 

. and also the surface-strand attraction alone is too weak for creating an adsorbed state, the combined 

' effect of the two attractions can lead to a naturated and adsorbed state. 

o : 

^ , PACS numbers: 82.35.Gh, 05.90.-l-m, 82.39.Pj 

in \ I. INTRODUCTION 

o . 

^ , The structure of DNA is the key for understanding its biological functioning, explaining why the physical features of 
' DNA have attracted attention over the last decades [1-7]. A known component of this structure is the Watson-Crick 
^ , double-strandedness: DNA is composed of two single-strand molecules lined up by relatively weak hydrogen bonds. 
The double-strand exists for physiological temperatures and is responsible for the stability of the genetic information 
stored in DNA. For higher temperatures the double-strand separates into two strands (denaturation). Many processes 
relevant for the functioning of DNA — such as transcription and replication of the genetic information and packing of 
^ ^ DNA into chromosomes — proceed via at least partial separation (denaturation) of the two strands due to breaking of 
hydrogen bonds [1, 2]. In addition, denaturation is important for a number of technological processes, such as DNA 
sequence determination and DNA mechanical nano-devices [2]. 

DNA denaturation is driven by changing the temperature or the solvent structure, e.g., the pH factor [1, 2]. There 
QQ ' are several generations of statistical physics models aiming to describe the physics of denaturation. Early models, based 
fSl , on the one-dimensional Ising model, focus on the statistics of hydrogen bonds modeling them as two-state variables 

■ (open or closed) [4]. More recent models describe a richer physics, in that they try to explore space configurations of 
Cn DNA [5, 5-9]. 

ly-j Most of the physics literature devoted to DNA denaturation studies this process in isolation from other relevant 

' processes involved in DNA functioning [3-9, 11]. However, denaturation is frequently only a component of a larger 
On , process, such as replication or compactification into a nucleosome, the basic structural unit underlying the chromosome. 

■ Here we want to study how another important aspect of DNA physics — adsorption of the double-strand DNA 
J> I on a surface — influences its denaturation. Surface adsorption of DNA is widely employed in biotechnologies for 

immobilization and patterning (drug or gene delivery) of DNA [10, 11]. There are in fact several pertinent situations, 
where both adsorption and denaturation of DNA are simultaneously at play. 
^ 1) For DNA at normal conditions (pH= 7 and NaCl concentration of 0.15 M) thermal denaturation occurs between 

■ ■ ' temperatures 67 C and 110 C (which are the temperatures for A-T and C-G unbinding, respectively) [1, 3, 4]. The 
denaturation temperature can be decreased by increasing the pH factor, i.e., by decreasing the concentration of free 
protons in the solvent, since the negatively charged phosphate groups on each strand are not screened anymore by 
protons and strongly repel each other. For the same reason, for the DNA adsorption on a positively charged surface, 
the increase of the pH will increase the electrostatic attraction to the surface. Thus at certain values of the pH factor 
and the surface charge, denaturation and adsorption may take place simultaneously. 

2) Surface adsorption can be realized by the hydrogen-bonding of the negatively charged phosphate residues to a 
negative surface (e.g., silica surface) [11, 13]. The effect is possible only when the electrostatic repulsion is sufhcienctly 
screened by the solvent cations. Thus the same factors (temperature, pH, solvent concentration) that decrease the 
inter-strand attraction, will weaken the DNA-surface binding [11]. 

3) The binding to hydrophobic surfaces (e.g., aldehyde-derivate glass, or micro-porous membrane) goes via partial 
denaturation which exposes the hydrophobic core of the double- helix and leads to the DNA-surface attraction [11]. 
Both naturation and adsorption are simultaneously weakened by increasing the pH [11, 12]. 
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4) Human DNA has a total length of 2 m bearing a total charge of 10* electron charge units. This long object is 
contained in the cell nucleus with diameter 10 /xm, which is comparable with the persistence length of DNA. Recall 
that the persistence length of a polymer is a characteristic length over which the polymer folds freely due to thermal 
fluctuations. For the double-strand DNA at normal conditions, the persistence length is relatively large and amounts 
to 50 nm or 100 base-pairs, while the persistence length of the single-stranded DNA is much smaller, about 1-2 nm 
(i.e., 2-4 base-pairs) [20]. This seems to create a paradoxical situation: not only the large, strongly charged DNA 
has to be packed in a very small compartment, but the DNA has to be replicated, repaired, and transcribed. The 
problem is solved by a hierarchical structure: the DNA double-helix is wrapped around positively charged histone 
(achieving partial charge neutralization), histones condense into nucleosome complex, which in its turn is contained 
in chromatin, etc. It was recently discovered that packing of DNA into nucleosomes with characteristic size much 
smaller than the persistence length of the DNA chain proceeds via transient denaturation of the double strand [21]. 
Denaturation reduces the persistence length and thereby facilitates the packing process. 

For all these processes we need to describe the DNA as a double-strand polymer interacting with an attracting 
surface. This will be the goal of the present paper. 

Needless to say that there is an obvious situation, where the double-stranded structure is not relevant for the 
adsorption. If the two strands are too tightly connected, their separate motion is not resolved. This case is well known 
in literature and — due to a large persistence length of a double-strand DNA — can be described via an effectively single 
semi-flexible chain interacting with the surface [18]. These studies complement the classic theory of the flexible chain 
adsorption, extensively treated in literature [3, 15, 16]. The electrostatic effects of the DNA adsorption, modeled 
via a single Gaussian chain, are studied in [22]. Another recent activity couples the Ising-Zimm-Bragg model for the 
helix-coil transition with the known theory of flexible chain adsorption on solid surfaces [17]. While interesting for 
their own sake, the results of Ref. [17] do not apply to DNA adsorption-denaturation, since the main assumption of 
Ref. [17] — that the helical pieces of the polymer interact with the surface much stronger than the coiled ones — does 
not hold for DNA. 

This paper is organized as follows. In section II we define the model we shall work with. It describes two flexible 
chains interacting with each other and with an attracting solid surface. Section II also recalls the known correspon- 
dence between the equilibrium physics of flexible polymers and quantum mechanics. In its flnal part this section 
discusses limitations of the studied model in including volume interactions. Section III recalls the quantum mechan- 
ical variational principle which will be the basic tool of our analysis. Section IV shows that if both polymers are 
adsorbed on the surface, they do not denaturate via a phase-transition. Section V discusses collective scenarios of 
binding, while section VI studies conditions under which the naturated and/or adsorbed state is certainly absent. 
The next section presents the phase diagram of the model. The last section summarizes our results. Some technical 
issues are discussed in appendices. The reader interested in the qualitative message of this work may study section II 
for learning relevant notations and then jump to section VII, which discusses general features of the phase diagram. 
A short account of the present work has appeared already in Ref. [14]. 



II. THE MODEL. 



When the motion of the single strands is resolved — i.e., when the inter-strand hydrogen bonds are relatively weak, 
as happens next to denaturation or unzipping transitions — DNA becomes a complex system with different, mutually 
balancing features at play. A realistic model of DNA should take into account the stacking energy between two base 
pairs and its dependence on the state (open or closed) of these pairs; helical structure of the double-strand; intra- 
strand and inter-strand volume interactions (e.g., self-avoidance); the pairing energy difference between A-T and C-G 
pairs (respectively, Sk^T and bk^T under normal conditions), etc. Such fully realistic models do not seem to exist; 
there are, however, various models with different degrees of sophistication which are intended to capture at least some 
features of the double-stranded structure [3, 5, 5-9, 15, 18]. 

The model we shall work with disregards almost all the above complex aspects and focuses on the most basic 
features of the problem. It consists of two homogeneous flexible chains interacting with each other and coupled to the 
surface described as an infinite rigid attractive wall. 

Consider 2A'^ coupled classical particles (monomers) with radius-vectors fi\i~ and f2\k {k = 1, ...,N) and potential 
energy 
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y {ra\k - ra\k-i) + V(f„|fc) 



(1) 



where — r^^, — f2\kj so that \ri\k — ^2\k\ is the distance between two monomers, U is the inter-strand potential and 
V is the surface-monomer potential. The harmonic interaction with stiffness K (Gaussian chain) between successive 



3 



monomers in each strand is responsible for the Hnear structure of the polymers. 

The system is embedded in an equilibrium thermal bath at temperature T = (/cb = 1)- The quadratic kinetic 
energy of monomers is irrelevant, since it factorizes from the partition function and does not influence the equilibrium 
probability distribution: 

(2) 

e-/3n(r,i.). (3) 

This model without adsorbing surface, i.e. V = 0, was mentioned in [15] and studied in [5] in the context of DNA 
denaturation. When the inter-strand interaction U{r) is absent, we get two independent flexible chains interacting 
with the solid surface, a well known model for adsorption-desorption phenomenon [3]. Recall that our purpose is in 
studying these two processes — i.e., surface-polymer interaction and inter-strand attraction — together. When taken 
separately, these processes are well studied and well understood. 

Note that for the considered Gaussian chain model the stiffness parameter K relates to the characteristic persistence 
length Ip as [3]: 

T 
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A. Specification of the surface-monomer potential. 

We assume that the surface can be represented as an infinite, solid plane at z = (the role of the solid surface can 
be played by any body of a smooth shape and the size much larger than the polymer length). Thus for the probability 
distribution (2) one has [for a = 1, 2 and k = 1, . . . , A''] 

P{ra\k) = P{Xa\k,ya\k, ^a|fc) = 0, for Za\k < 0- (5) 

This boundary condition should be imposed as a constraint in (2). 

The remaining part of the surface-strand interaction is described by a negative (attractive) potential V that depends 
only on the third coordinate: V{fa\k) = V{za\k)- The potential V{z) will be assumed to be short-ranged: it is negative 
for 2; ^ and tends to zero sufficiently quickly for 2; ^ 00. 

Let us continue the specification of the potential V(z) taking as an example the electrostatic attraction between one 
negatively charged DNA strand and a positively charged surface; see, e.g., [16]. We denote by a the surface charge 
density, q stands for the monomer charge (for DNA the effective monomer charge is roughly q le, where e is the 
electron charge) and e is the dielectric constant of the medium into which the polymer is embedded (e « 80 for water 
at room temperature). Now the interaction energy between the surface area dxdy and one monomer reads: 

— e-'^^Majdy, (6) 
er 

where r = ^/x^ +y'^ + z"^ is the distance between the surface area dxdy and the monomer, while fco is the inverse 
Debye screening length. This interaction leads to attraction for opposite charges: qa < 0. The full expression of the 
inverse Debye screening length is well-known: 



kn = 2nlBj2^aZl ^B = ^, (7) 
y a 

where Ib is the Bjcrruin length, and where Ua and Za are, respectively, the concentration and valency of ions of the 
sort a present in the solvent (so that the ion charge is Zac). The summation in (7) is taken over all sorts of ions 
present ^ . Under normal conditions the Bjerrum length is ~ Inm. At this length the electrostatic interaction energy 



^ The quantity ^ X^u^a-^a called ionic strength. 
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becomes comparable with the thermal energy T. The Debye length l/ku varies between ~ 0.5nm and ^ 1.5nm under 
physiological conditions. For pure water it is much longer: l/fc^ ~ lfj,m. 

Integrating (6) over x and y from — oo to +00, we get for the surface-monomer interaction [16]: 

V(.) = ^e-^-. (8) 
eku 

Thus the strength of the potential is while the inverse characteristic length is (expectedly) 1/ku- The potential 

V{z) is short-ranged for all other relevant mechanisms of adsorption (hydrogen-binding, hydrophobic interactions, 
cation exchange). This means, in particular, that dzV{z) is finite for all these mechanisms [23]. 

Returning to (8) we note that for a single flexible polymer interacting with the surface the adsorption problem was 
solved in Ref. [16] within the Schrodinger equation approach to be discussed below in detail; see in particular Eq. (23). 
The adsorption-desorption phase transition transition temperature found in [16] reads: 

8.33 7r|(T(j'| 

where Ip is the persistence length from (4). 

Let us estimate the Debye length as k^^ ~ Inm, the single-polymer persistence length as Ip ~ Inm, and assume 
that the surface contains Z elementary (electron) charges per 1 nm x 1 nm. Normally Z ^ 1, though stronlgy charged 
surfaces achieve Z = 10 — 20. Taking the effective monomer charge one elementary charge (which is a typical value for 
a single-strand DNA) and recovering the Boltzmann constant, we see that (9) predicts Tc of order of room temperature 
(300 K). 

When looking at concrete parameters in (9) we should also recall that Eqs. (6-8) account for the surface-monomer 
electrostatic interaction, while the monomer-monomer electrostatic interaction within the single polymer is neglected. 
This is possible when the surface charge al^ at the area Ip x Ip (where Ip is the persistence length of the single strand) 
is larger than the monomer charge: 

(10) 

This condition will be satisfied for strongly charged surfaces Z « 10. 



B. Specification of the monomer- monomer interaction between the two strands. 

The inter-strand potential ZY(|rin, — r2|fe|) collects the effects of hydrogen-bonding, (partially) stacking, and possible 
electrostatic repulsion. We again assume that it is purely attractive, short-ranged and goes to zero sufficiently fast 
whenever the inter-particle distance — r2|fe| goes to infinity. In particular, the short-ranged features implies that 

drr^U{r) is finite. 

Several concrete examples of the inter-strand potential U were studied and favorably compared with denaturation 
experiments in [5-7]. For example, Ref. [6] studies the Morse potential 

l({r) = ue-^'^ie-"'' - 2), (11) 

where 1/ is the potential strength and ^ is its characteristic range. Within the Schrodinger equation approach [see 
(23) below] Eq. (11) predicts a second-order denaturation transition at the critical temperature: 



(12) 



Note that the appearance of the factor a^l^ in (12) is similar to the that of the factor k^l^ in (9). Here are the 
standard estimates for the parameters in (12): z/ ~ O.OleV, Ip ~ Inm and alp ~2 [7]. These produce from (12) 
Tc ~ 400K, which by the order of magnitude coincides with experimental values [7]. 



C. Effective Schrodinger equation. 



It is known, see e.g. [3, 7], that in the thermodynamical limit TV ^ 1 the free energy of flexible polymer in an 
external potential is determined from an effective Schrodinger equation; see Appendix A for more details. A sufficient 
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condition for validity of the Schrodinger equation approach is that the characteristic length D over which the polymer 
density changes is much larger than the persistence length Ip-. 

D » Ip. (13) 

This condition is always satisfied in the vicinity of a second-order phase-transition, where D is naturally large for a 
fixed Ip. If condition (13) is satisfied for a short-range potential — see (8, 12) for relevant examples — this potential is 
necessarily small for those distances ~ D, where the flexible polymer is predominantly located [3]. 
For the considered two-strand situation the Schrodinger equation reads 

= E^-, (14) 

^ = + ^(-^«)] + (15) 

where [using also (4)] 

V{z) = Kp^V{z) = ^V{z), U{r) = K(3^U{r) = ^U{r). (16) 

ip ip 

If there is a gap between the lowest two eigenvalues of H, the ground state wave-function 4" determines the monomer 
statistics as 

n(rl,r1) = *2(ri,f2), (17) 

where n(ri, r^) is the probability distribution for two neighboring monomers on the strands for the considered trans- 
lationally invariant system. 

Recalling the known correspondence between the flexible polymer physics and (stationary) quantum mechanics [3] , 
we can think of ri,2 = {xi,2, 2/i,2, ^1,2) as the position vectors of two quantum particles representing the strands, while 
r = r*! — r2 is their mutual position. 

The eigenvalue E is the energy of the quantum pair. It is related to the free energy fN of the system as 

E = 0Hf + 3plln—. (18) 

Since the surface is described by an infinite potential wall, we have the following boundary condition for the wave 
function ^ 

*(ri,f2) = 0, if zi<0, or ^2 < 0. (19) 

Both V{z) and U{r) are attractive, V < 0, U < 0, and short-ranged, that is dzV{z) and drr^U{r) arc finite. 
When [/ = 0, the Hamiltonian H reduces to two uncoupled strands (or two uncoupled quantum particles), each one 
in the potential V{z). The corresponding Schrodinger equation for the ^-coordinate of one strand reads from (15): 

[- ^d^ + Vizmz) = Ei^iz), i;{z = 0) = 0. (20) 

It is well-known that if V{z) is shallow enough, no bound (negative energy) state exists, while the second-order 
binding transitions corresponds to adsorption of a single flexible polymer [3] . The physical order-parameter for this 
transition is the inverse square average distance from the surface, l/(z^), which is finite (zero) in the adsorbed 
(desorbed) state. It is useful to denote by /x the dimensionless coupling constant of ^ = iiV such that (for [7 = 0) 
the adsorption threshold is 

Mc.o = 1- (21) 



^ In fajct, one should be more careful, when defining the boundary condition (19). For the two-particle case it appears to be necessary to 
fix not only the continuity of and its value at the wall, as Eq. (19) docs, but also the behavior next to the wall: one has to require 

that when zi and Z2 go to zero simultaneously, ^ oc ziZ2. Otherwise, there will be (continuous) wave- functions which provide a bound 
state for two-particles with an arbitrarily weak y < and arbitrary weak inter-particle interaction U < 0, though the single particle 
needs a critical strength of V to get into a bound state. This obviously pathological situation is prevented by the additional boundary 
condition ^ oc 2122. For the wave- functions we shall consider below this additional boundary condition will be satisfied automatically. 
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Note that the adsorption of a single strand DNA is a part of the renaturation via hybridization [2] , a known method 
of genetic systematics. 

For the example (8) the concrete expression for fj. reads from (9): 

Analogously, switching off both V{z) and the wall, we shall get a three-dimensional central-symmetric motion in 

the potential U (r) which again is not bound if U is shallow. This second-order unbinding transition with the order 
parameter l/(r^), where r is the inter-strand distance, corresponds to thermal denaturation (strand separation) of 
the double-strand polymer [1, 2, 5, 7]. 

The Schrodinger equation for the radial motion in the absence of the surface reads from (15) [7] 

[~d^ + Uir)]xir)=Exir), x{r = 0) = 0, (23) 
where %('') is related to the original wave- function as 

V)(r) = (24) 

Note that (23) is again a one-dimensional Schrodinger equation, but as compared to the equation (20), Eq. (23) 
contains an additional factor ^ next to the kinetic-energy term d^. This factor arises due to effective mass; see [23] 
for more details. 

Let us write likewise U = AC/, where A is the dimensionless naturation strength. We take the naturation threshold 
in the bulk to be 



For the example (11), A reads from (12): 



Ac,o = 1. (25) 



When the wall is included, i.e., condition (19) is imposed, the strands loose in the adsorbed phase part of their 
entropy. This is known to lead to a fluctuation induced effective repulsion [19]. 

Let us now remind that the physics of weakly bound quantum particles does not depend on details of binding 
potential [23] . Thus for qualitative understanding of the situation one may employ the delta-shell potential, which is 
easily and exactly solvable and has very transparent physical features; see Appendix B. 

D. Relevant coordinates. 

Let us now return to the basic equation (20). It is convenient to recast this equation in new coordinates: 

vi = ^{xi + X2), V2 = ^{yi + 2/2), (27) 
xi - X2 = pcosip, yi - y2 = p sirup, (28) 

where 

< p, < ^ < 27r, (29) 
and to re-write the Schrodinger equation (14, 15) as 



1 (2 d d Id 1 1 52 , TT./ X 



Z2) + U{\fi-f2\)}^ = E^. (30) 
It is seen from (30) that the variables separate, since ^'(ri, r2) can be written as 

*(ri,f2) =tp{p,Zi,Z2)^l{vi) ■ip2{v2)'tp3{>P), (31) 
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and the lowest energy levels is to be found via the following equation 

+ 5 + 5} * + {"^^"^ + ^^"^ + f/(y7T(i7^)} ^ = (32) 

Thus due to the translational invariance along the surface and the invariance under rotations around the z-axis, 
we are left with three independent coordinates: the projection p of the inter-particle distance on the surface, and the 
distances zi , Z2 between the particles and the surface. 

Note that within the quantum mechanical setting the problem described by (32) corresponds to a three-body 
problem, where the role of the third body (with infinite mass) is played by the surface. 

E. Common action of the surface-strand and inter-strand potentials. 

In this and subsequent subsection we shall discuss two possible limitations of the present model. 

Above we combined together the surface-strand interaction potential V, which was derived separately from studying 
interaction of the surface with one flexible strand, and inter-strand potential U deduced from studying two flexible 
strands without the surface. While this type of combining is widely applied in all areas of statistical physics, its 
applicability needs careful discussions in each concrete case. For instance, it is possible that the presence of adsorbing 
surface will directly influence the inter-strand potential. Let us discuss one (perhaps the major) example of that type 
pertinent for the studied model. 

It is well-known that the two strands of DNA are negatively charged [1] . For the double-stranded DNA under normal 
conditions the inter-strand repulsion is screened by positive coimterions, so that the hydrogen bonding can overcome 
the electrostatic repulsion and create an effective attraction, which is then the main reason of inter-strand binding 
[1]. Once DNA denaturates and separates into two strands, the counterions are released into the ambient medium 
and are clouded around each strand. However, for temperatures not very far from the denaturation temperature the 
counterions continue to screen the electrostatic repulsion, so that once the temperature lowers below the denaturation 
transition temperature, the two strands reversibly assemble back into the double-strand [1]. We stress that the fact 
that (partially released) counterions still provide a sufficient screening follows from the existence of the observed 
reversible renaturation transition. 

When DNA denaturates in the presence of a positively charged siirface the cloud of screening counterions around 
each strand will tend to rarefy. This will increase the screening length and make the overall inter-strand interaction 
repulsive. However, this is possible only for strongly adsorbed strands, where the majority of counterions are within the 
direct inffuence of the surface charge. In the present work we focus on weakly bound strands, where the characteristic 
length of the adsorbed layer D is much larger than the persistence length Ip (approximately Inm in normal conditions), 
which is of the same order of magnitude as the Debye screening length l/kn; see (13). Thus the majority of counterions 
will not feel the adsorbing surface, and in this case we do not need to account directly for the influence of the surface 
on the inter-strand potential. For strongly adsorbed DNA strands, i.e., for D ~ I/^d, it can be necessary to couple 
directly the inter-strand potential with the degree of adsorption. 

F. Self-avoidance and of electrostatic volume interactions. 

In Hamiltonian (1) we accounted for the surface-strand and inter-strand interaction, but neglected all the volume 
interactions such as self-avoidance and (for charged polymers) electrostatic interaction between varioiis monomers. 
It is important to note that the volume interactions coming from the intra-chain contributions can be accounted for 
within the present model via renormalizing the persistence length Ip-, see (1) and (4) for definitions. As shown in 
[25] for a single flexible polymer interacting with electrostatically adsorbing surface, the self-avoiding interactions and 
electrostatic volume interactions renormalize the persistence length. Provided that the Debye screening length 1/kj) 
is not very large — a sufficient condition for this is kr>lpVN 1, where N is the number of monomers [25] — both 
self-avoiding and electrostatic volume interactions lead to an effective persistence length Ip, which differs from the bare 
persistence length mainly by the factor A^^/^°: Ip ~ 7V^/^° Ip [25] (the remaining part of renormalization is numerical 
factors, which are not essential for the present qualitative discussion) ^. Once the persistence length is renormalized, 



^ In more detail, Ref. [25] considers a continuous polymer model with length L and reports for the square of the effective persistence 
length Ip ~ L^/s Pqj. ^^e present dsicrete model we take naturally L oc TV. 
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one can still use the flexible polymer coupled to an adsorbing surface [25] . Thus the transition temperatures (9) and 
(12) are divided by factor N^/^ , where N is the number of monomers. Now for the typical single-strand DNA length 

~ 10^ this renormalization will not make any substantial change in transition temperatures, though it is essential 
for longer polymers, N > 10^. In particular, for such a long polymer the persistence length may increase to an extent 
that the condition (13) will be violated. 

We will see below that for qualitative conclusions of this paper, the precise form of the renormalized persistence 
length is not essential, provided that one can still employ the Schrodinger equation (20) for describing denaturation and 
desorption. The main reason for this is that the renormalization of the persistence length homogeneously renormalizes 
both dimensionless couplings A and /i in (25) and (21), respectively. 

The above discussion does not account for the inter-chain volume interactions and thus should not create an im- 
pression that the full volume interactions effect for two coupled chains can be described via a renormalized persistence 
length. It is clear that one needs a more specific study of volume interactions for the present model. Since such a 
study poses immence analytical problems, it will be concluded at a later time. 



III. VARIATIONAL PRINCIPLE AND THE EXISTENCE OF THE OVERALL BOUND STATES. 



Note that Eqs. (32, 19) follow from a variational principle: 

5J{^} = 0, 



(33) 



with 



nOO nOO nOO 

^{V'} = / / / pdpdzidz2 

Jo Jo Jo 



(34) 



where ijj is taken real, since we are interested in bounded (discrete-level) states. We already assumed that ip is properly 
normalized: 



JO 



pdpdzidz2-tp^izi,Z2,p) = 1. 



(35) 



If either V{z) = or [/ = 0, the criterion for the existence of a bound state is well known, since it reduces to the 
existence of a negative energy in the spectrum, or equivalently to the existence of a physically admissible (satisfying 
the proper boundary conditions) wave-function with a negative average energy. 

The situation is slightly more delicate when the two potentials V and U act together. Let us assume that either 
for V{z) ^ or for U {\fi — f2\) there are negative energy states. Denote by 



E{U}<0, E{V{zi) + V{z2)} = 2E{V} <0, 

respectively, the corresponding lowest (most negative) energies. 
Then it suffices to have a normalized wave-function tp with 

J{^} < E{U} + 2E{V}, 

for at least one overall bound, i.e., adsorbed and naturated, state to exist. 



(36) 



(37) 



IV. ABSENCE OF DENATURATION PHASE-TRANSITION FOR ADSORBED STRANDS. 

Let us return to the variational principle (34) and assume that V{z) is strong enough to create at least a single 
(lowest) bound state with energy E{V} < 0. Denote by (f){z) the corresponding lowest-energy normalized wave 
function: 

-lct>"{z) + V{zmz) = E{V} ct>{z). (38) 
For the overall problem we shall employ the following variational wave-function: 

-tpiP, Zi,Z2) = 4){zi) (t){z2) iip), (39) 
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where ^{p) is an unknown, tentatively normalized, viz. 

3 

dppeip) = 1, (40) 



/ 

Jo 



wave- function, to be determined from the optimization of (34). Note that in (39) the boundary conditions for the 
surface are satisfied via ^{zi) (p{z2)- 

Substituting (39) into (34) and varying it over ^, we get an effective Schrodinger equation for ^(p): 

--^P-^}^ + {UM~e}i = 0, (41) 



p dp dp 

where Uee{p) is an effective potential: 

POO POO 

Ues{p)= dzi dz2 (p^zi) <j)^{z2) U{ Vp' + {zi - Z2)^ ), (42) 
Jo Jo 

and where e is the reduced energy 

e = E- 2E{V}. (43) 

Two main point about the effective potential (42) is that it is attractive (since so is U) and goes to zero for p oo. 
The last feature follows from the analogous one of U (r) and the fact that 4i{z) are normalizable. A more explicit form 
for Uqs can be obtained by assuming that U{r) is a delta-shell potential 

U{r) = --d{r-ro), (44) 
ro 

with the strength A > and the attraction radius ro > 0. The transparent properties of this potential are recalled in 
Appendix B. The critical binding strength of this potential is 

Ac,o = 1, (45) 

as given by (B9). [When comparing Eq. (44) with Eq. (B2), note that the additional factor 2 comes from the reduced 

mass.] 

Using (44) we now obtain from (42) after changing variables: 



Uesip) = - 




(46) 

It is now seen explicitly that Ues{p) is zero for sufficiently large p. 

Note that (41) has the form of two-dimensional Schrodinger equation for an effective particle in the attractive 
potential Ues{p). It is well known that any (however weak) attractive potential in two dimensions creates a bound 
state [23]. Thus there is a normalizable function ^{p) such that e in (42) is negative. This means that 

E < 2E{V}, (47) 

and, according to our discussion in section III, there is an overall bound (naturated and adsorbed) state provided 
V{z) creates a bound state. In our model a sufficiently attractive surface potential confines fluctuations of the two 
strands and prevents the dcnaturation phase-transition (this however does not mean that the denaturation is absent 
as a physical process; see below). 

The physical reason for the existence of an overall bound state for an arbitrary small potential is a peculiar two- 
dimensional effect: the weakly singular attractive oc l/p^ potential [19] ^. Indeed changing in (41) the variables 



* One-dimension in this respect is not much different from the three-dimensional situation. The known statement on the existence of 
bound state for any small one-dimensional potential is connected with a different mechanism, that is, with allowing all values of the 
one-dimensional coordinate (no infinite wall at the origin). The two-dimensional situation is indeed peculiar in this respect. 
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as 



we get 



^ + {uM-^-s}^ = 0. (49) 



Eq. (48) implies 



m = 0, (50) 



i.e., the existence of the infinite wall at p = for the effectively one-dimensional Eq. (49). ft is, however, seen from 
(49) that there is also an attractive potential l/(4p^). ft is known that if the strength of such a potential is larger 
than 1/4 seen in (49), the (effective) quantum particle will fall to zero, i.e., the ground state will be minus infinity [23]. 
The value 1/4 is just at the border of this phenomenon and, therefore, any attractive short-range potential acting in 
addition to l/(4p^) suffices to create a bound state [19]. 

To illustrate the behavior of UeS for weakly bound state of the potential V{z), let us assume that V{z) is also a 
delta-shell potential: 

V{z) = -£-5{z-zo). (51) 

Recall that we still have an infinite wall at 2; = and that for the delta-shell potential the bound state exists for 

M > Mco = 1, (52) 

see Appendix B for details. If now p. is close to one, the energy E{V} = —k'^/2 is small. Working out (46) with help 
of Eq. (Bll), which essentially reduces to 



(/)(z) oc \/2fc e-'=^ (53) 

we get 

t/eff(p) = -2A^%=4 (54) 

V'^o - P 

Since for small k, the wave- function ip{x) is almost delocalized, the efii'ective potential Ues{p) is proportional to k and 
goes to zero for fc — > that is for z/ — > 1. In other words, the trial function (39) does not predict any (overall) binding 
for 

M < 1- (55) 

Note however that although for > 1 any inter-strand attraction is able to prevent the denaturation phase 
transition, the energy s in (49) is exponentially small for small Ugs, i.e., small A or small k. Recall that this energy 
is estimated as [23] 



2 

e -5- exp 

^0 



dppUes{p) 



4e-V(^'='-«). (56) 



r, 







Thus for a small A or A; we get a very large separation between the strands. In this sense the (incomplete) 
denaturation phenomenon without the phase transition is present in our model. 

In summary, the main physical message of this section is that if the two strands are localized near the surface, the 
overall DNA molecule does not melt via a phase-transition with increasing the temperature: there is only a smooth 
crossover from tightly bound to a (very) weakly bound state. The cause of this effect is that the surface confines 
fluctuations of each strand. Mathematically this is expressed by an additional attractive potential —4^ in (49). 

This result was obtained without taking into account various realistic features of DNA. It is possible that the 
denaturation transition in the adsorbed phase will recover upon taking into account some of those neglected features, 
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e.g., volume interactions between the two strands and within each strand (see [9] for a prediction of such a transition 
in a different model of DNA that partially accounts for volume interactions). 

We nevertheless expect that the obtained result will apply, at least qualitatively, to denaturation-renaturation 
experiments, and will be displayed by facilitation of the naturation in the adsorbed phase. We are not aware of 
any specific experiment done to check the renaturation-facilitating effect of an attractive surface. There are, however, 
somewhat rcilatcid experiments showing that the renaturation rate can significantly increase in the condensed (globular) 
phase of single-strand DNA [24] . This condensed phase is created by volume (monomer-monomer) interactions. The 
effect was obtained under rather diverse set of conditions, but to our knowledge it did not get any unifying explanation. 
The analogy with our finding is that in the condensed phase fluctuations of the single strand DNA are also greatly 
reduced as compared to coil (free) state. 



V. COLLECTIVE BINDING. 



With the aim to understand the situation when V{z) alone does not provide any binding, we take for the variational 
function 

i^{zi,Z2,p) = <j>{zi,Z2)£.{p). (57) 

As compared to (39) we do not require that zi and Z2 are factorized, and we are going to optimize over (jji^ii ^2)- In 
contrast, ^(p) is a fixed, normalized (see (40) ) known function. 
Substituting (57) into (34) and varying over (j){z\,Z2) we get: 

"Hfi'^fi}'^^ ^^^^'^ ^ ^^^'^ ^ ^"^^^'^ ~ ^''^ - i^i} <^ = 0, (58) 

where 

V,s{z)= / dppe{p)U{^/^T^). (59) 
Jo 

Recall that by the very meaning of the variational approach Ei provides - -for any A and any normalized function 
^(p) — an upper bound for the real ground state energy. Eq. (58) describes two one-dimensional particles with 
inter-particle interaction ycfi(l-Zi — Z2I) and coupled to an external field V{z). 

For the inter-particle interaction given as in (44), this effective potential Ven{z) reads 



V,s{z) = -Xe{ro ~ z)e (^^Jrl-z^^ 



(60) 



We are now going to show that Eq. (58) predicts binding — that is, it predicts Ei < and a localized normalizable 
wave-function (j){zi, Z2) — at the critical point p, = 1 of the potential V{z). To this end let us calculate the perturbative 
correction AE introduced by the effective potential Vca- At first glance the application of perturbation theory is 
problematic, because we search for a nearly degenerate energy level. However, due to strong delocalization of the 
corresponding wave-function, the matrix elements of the perturbing potential Ves appear to be small as well, and 
applying perturbation theory is legitimate. This will be also underlined below by a perfectly finite behavior of the 
second order perturbation theory result. 

Recall that in the first two orders of the perturbation theory we have [23] 

AE^E,- 2E{V{z)} = (OiKff |0) - j^^ dK f^^f^^^^y (61) 
2E{V} = -k^, (62) 

where (zi,Z2|0) = (j){zi)(j)(z2) is the lowest energy state of the unperturbed system, and where the integration over 
K involves all excited wave-functions of the unperturbed two-particle system with wave- vector K and energy ek (all 
these wave- functions are in the continuous spectrum). Note that there are three orthogonal families of these states: 

(t>{zi)(l){nz2,n), ^"^"2 2"' 

(j){z2)(t){nzi,n), = ^ 2"' ^^^^ 

2 1 

(j){nzi,ni) (t>{nz2,n2), e„i„2 = + (65) 
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where (j){nz, n) are the corresponding single-particle excited (continuous spectrum) wave-function with the wave- 
number n. These wave-functions are normalized over the wave-number scale; see Eq. (B14) in Appendix B. This type 
of normalization is important for the integration over the wave- number K in (61). 

The first-order contribution to Ai? appears to be zero for A; ^ 0+ (i.e., for /x 1+). Indeed, we can use (53) for 



ci>{zuZ2)=(l>{zi)cl>{z2) = 2ke-''^''+'^\ 



(66) 



to conclude 



OO /'OO 



(0|i/|0)= / / dzidZ2V;ff(|zi-Z2|)(/.2(zi)02(z2) 
JQ Jo 



POO /'V 

I dv / duVes{u)i 
Jo 



V + u 



2 y ' V 2 

Using (63-66) we shall calculate various matrix elements entering into (61): 



2k dwe-2'='' / duV^s{u) = 0{k). (67) 
Jo Jo 



OO /'OO 



{0\H\n) = / dzidz2(l)izi)(l){z2)Vef[{\zi - Z2\)(t>izi)(l){nz2,n) 
Jo Jo 

' n{u + v) 



/>00 f^V 

V2k2k dv due-'''^'^^+^^^^Ves{u)(j) 
Jo Jo 



(68) 



{0\H\nu ns) = 

= 2k dv £ du e--/^ VMu) (^^^i^ 
This results to the following formula for AE, 



dzi dz2(j){zi)(f>{z2)Ves{\zi - Z2\) <f>{nzi,n) (f>{nz2,n) 

n2{v - 



ni 



,n2 



AE = -2 dn '^j ' - / dnidn2 'J ''^ . 
Jo ^ + V Jo !^ + !^ + ^ 



(69) 



(70) 



Working this out and going to the limit A; ^ 0+ (i.e. 



AE = -i 



+ 



f 

Jo 



duVcs{u) 
dni dn2 



2 + ni + n^ 



1"*") we obtain 
dn 



^~a„->(-„)^(- „)]%<„. 



(71) 

(72) 



This expression for AE is finite in the limit fc — > (see Appendix C for details), and proportional to the squared 

perturbation strength [Jp°° du T4fl (?i)] . In the limit ji ^ 1 (and for sufficiently small A) wc arc in the situation where 
neither V{zi) + V{z2) nor U alone create bound states. Recalling our discussion in section III on the existence of 
bound states as reflected in the magnitude of variational energy, we conclude from AE < that the present approach 
does predict binding for = 1 and for sufficiently small A. Since the ground state is supposed to be continuous, the 
very fact of having a negative energy for fJ- = I and not very large A implies that a bound state will exist for 



Ac > A > 0, 



/Xc < M < 1; 



(73) 



where neither of the potentials V and U alone allows binding. Here Ac > 1 is some critical value at which the real 
ground-state energy is equal to E{U}; recall our discussion in section III. Note that the precise form of ^(p) is 
irrelevant for the argument. This function has to be normalized and such that the effective potential VeS does not 
become large for a sufficiently small A (and, of course, does not vanish for a finite A). For the rest it can be arbitrary. 

Thus in view of (73) we have found an example of so called Borromean binding, where the involved potentials do 
not produce bound states separately, but their cumulative effects lead to such a state. It is seen from (57) that this 
unusual type of binding is connected with correlations between the z-components of each particles and separately 
with correlations between their x and y components (which enter via p). 

Note that for three (or more) interacting point- like particles (instead of two particles and a surface) this effect was 
predicted in nuclear physics; see, e.g., Ref. [26] for a review. 
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VI. NO-BINDING CONDITIONS 
A. First method. 

Here we shall consider certain lower bounds on the sought ground state energy. Although these bounds are basically 
algebraic, they are non-trivial, and they allow to find out under which conditions both the adsorption and naturation 
arc absent . In this way wc complement the study of the; previous section. We employ — with necessary modifications 
and elaborations for our situation — the method suggested in [27]. 

Note from (19) that the presence of the infinite wall can be modeled via the boundary condition at the plane z = 0: 

*(n,f2) = 0, if zi=0, or 2:2 = 0. (74) 

Though the physical content of the problem demands that 5'(fi,r2) is also zero for zi <0 ov Z2 < 0, we can formally 
require only (74) and continue the potential V{z) to z < via 

Vi-z)^Vi\z\). (75) 

The ground state energy of the new problem defined with help of (74, 75) will be obviously equal to the ground state 
of the original problem. 

Let us now introduce a Active particle with the mass M and the radius vector 

r3 = {X3,y3,Z3)- (76) 

Now Eq. (75) is generalized to the corresponding translation-invariant interaction with the fictive particle: 

y(|zfe-Z3|), fc = l,2. (77) 

It is again obvious that upon taking the limit M ^ 00, the motion of the fictive particle will completely freeze, rs 
will reduce to a constant which can be taken equal to zero. 

Thus the three-particle (two real particles plus the fictive one) Schrodinger equation reads analogously to (14, 15) 

+ V(|zi - Z3I) + Vi\z2 -Z3\) + U{\ri - ral) - E{M) \^ = 0, (78) 



2M df^^ 2 dfi^ 2 9r2^ ^' ^' "''^ ^' " ^ ' } 

the correct two-particle energy being recovered in the limit M — » 00. 
Note that the boundary conditions (74) are modified as well 

= 0, if zi = Z3, or Z2 = Z3. (79) 

It is seen that the Hamiltonian in (78) is invariant with respect to simultaneous shift of all three radius vectors 

Tfe (fc = 1,2,3) by some vector. Since we consider a finite-particle quantum system, symmetry of the Hamiltonian 
implies the symmetry of the corresponding ground-state wave-function. Thus we deduce for this function 

* = «'(rl - r2, rl - rg, ra - rj), (80) 

which implies 

Wc shall now decompose the Hamiltonian in (78) such that (81) is employed and that the separate sectors of the 
problem — i.e., surface-particle and intcr-particlc interaction — are made transparent: 



H = -77i7 ^-oT7^-oT7^ + ^(1^1 - ^3|) + V{\Z2 - Z3\) + C/(|r-i - rl,]), (82) 



1 1 92 I 

2M dfs 2 dfi 2 df2 
= Hq + H\2 + -ffi3 + H23 , (83) 

H =--(— — —\(a— b— b—^ (84) 

° 2 \ 9r3 dfl df2 ) \ dfi, dfl dfi J 



c 



d X d 



2 



^---HrT^^-TT^a;^' +ni^i-3i) (85) 

if23 = -^fr^^-T^,^) +V{\Z2-Z^\) (86) 
^-^-^'^ (29^-29^1 +^(l^-^^2|). (87) 



2 \ 1 -h a; df2 1 + x dr; 
Id 19"^ 
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The coefBcients a, b, c and d are read off directly from (82-87): 



« = -7^^^. (88) 



2^ 

(1 + 2xy 

(l + x)^ 
(1 + 2x) 



2 • 



(89) 



2x(l+x) , , 

where the hmit M — > oo has already been taken. Here a; is a free parameter; the boundaries of its change are to be 
determined below. 

Let us now take average of the Hamiltonian H with the ground-state wave function The term {^\Hq\^) is 
zero due to (81). Wc shall now establish when the remaining terms in are certainly positive, that is when 

bound — i.e., naturated or adsorbed — states are certainly absent. 

Changing the variables as 



a = (1 + 2x)rx + fa, 6 = xr^ + ^3, (91) 



one reduces i?i3 to a form 



= -f7z^ + ^(16^-26.1), (92) 

where and arc the third components of the vectors ^1 and ^3, respectively. The constant factor 2^3^ will 
obviously not change the binding conditions. Recalling boundary conditions (79) we see that i^\H\s\^) is certainly 
positive for 

/X < c, (93) 

where \i is the coupling constant of V , such that Hy^ with c = 1 has the binding threshold /x = 1 [compare with (21, 
52)]. Obviously, (^|if23|^) is positive under the same condition (93). 
As for (^|iJi2|^) we change the variables as 

r\i = ri - r2, R\i = — - — , (94) 



to see that H12 takes the form 



Thus, (^'lifisl^) is certainly positive for 

A < 2d, (96) 

where A is the coupling constant of U, such that H12 with 2rf = 1 has the critical binding threshold A = 1 [compare 
with (25,45)]. 

Let us now recall that we employed c and 2d as inverse effective masses which should be positive; thus we should 
restrict ourselves to the situations 

X > (97) 

and X < —1, as seen from (89, 90). As inspection shows, the relevant no-binding condition is produced for x changing 
from zero to plus infinity, i.e., for the branch (97). 

Thus, under conditions (93, 96), where the limit M ^ 00 is being taken, the overall bound states are certainly 
absent. 
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B. Second method. 



Let us now turn to another, simpler way of deriving no-binding regions. For some range of parameters the present 

method will have a priority over the considered one, and then by combining the two methods we shall get an extended 
no-binding region. We return to the very original quantum Hamiltonian in (15) and write it as 

a f d d 
4 \ dfi df2 
a f d d 

1-a 



2 arl^ 
l-g d"^ 
2 df*i 



(98) 

2 

+ C/(|ri-f2|) (99) 
Vifx) (100) 
y{ri), (101) 



where 



< a < 1. (102) 

The term in (98) is seen to be always positive; for the term (99) we change the variables as in (94), to get that it 
is always positive for 

A < a, (103) 

while the terms in (100, 101) are both positive under 

< 1 - a. (104) 

Here A and \i have the same meaning as in (93, 96). Thus no binding is possible if (103) and (104) are satisfied 
simultaneously. 



C. Convexity argument. 

As for the last ingredient of our construction, we note that the coupling constants /i and A enter into Hamiltonian 
ff(/ti. A) in the linear way, and that the following convexity feature is valid for the ground state as a function of \i and 
A: 

min \H{viix + (1 - z^)M2, v\x + (1 - v)\i) ] (105) 
= min [!/if(/ii,Ai) + (l-z/)iJ(/i2,A2)] > i^min [iJ(/ii, Ai)] + (1 - z/) min [if(/x2, A2)] . (106) 

In other words, if in the phase diagram the binding — i.e., naturation or adsorption — is prohibited at points (/xi, Ai) 

and (/i2,A2) — that is min [ff(/xi,Ai)] > and min [ff(/X2,A2)] > — then there is no binding on the whole line 
connecting those two points, because from (105, 106) one has min \H.(y\x\ + (1 — v)\i2i v\\ + (1 — v)\i) ] > 0. 

Thus we draw together the bounds (93, 96, 103, 104) — under conditions (97, 102) determining the ranges of the 
parameters x and a, respectively — and complete it to a convex figure ensuring that for every two points belonging 
to (93, 96, 103, 104) the line joining them is also considered as binding-prohibited. The result is presented in Fig. 1. 
It is seen that there is the critical strength /ic = 0.25 — which is necessary for binding. 

The latter value of is special for the following reason: for A ^ 00, i.e., when the inter-particle attraction is too 
strong, the two particles are tightly connected to each other. The mass of the composite particle is two times larger, 
and (at the same time) the potential acting on it is two times larger. This leads to the adsorption threshold 11 = 0.25, 
which is independently obtained via the above no-binding conditions. 



VII. PHASE DIAGRAM. 



We are now prepared to present in Fig. 1 the qualitative phase diagram of the model. The axes of the phase 
diagram are A and /x. The dimensionless parameter A enters into the inter-strand interaction energy XU, such that 
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FIG. 1: Schematic phase diagram for the inter-strand coupling A versus the strand-surface coupling /i. The bold lines confine 
three thcrmodynamical phases. ND: Naturation and desorption. NA: Naturation and adsorption. DD: Desorption and 
denaturation. The critical naturation strength in the bulk is Ac = 1, for single strand adsorption it is Hc = 1. The following 
subregions are confined by normal lines, a: Domain described by the no-binding condition of section VI. b (bounded by the 
bold DD-NA lino and two straight segments): Borromean naturation and adsorption, c: Adsorption and naturation due to 
overcritical coupling to the surface. 



without the adsorbing surface the naturated phase of the two strands exists only for A > 1; see sections II B and II C 
for details. In this phase the two strands arc localized next to each other and their fluctuations are correlated. The 
typical form of A for the considered short-range potentials is 

^=^' « 

where c is a numerical prefactor, is the strength of the inter-strand potential (i.e., the modulus of its minimal value), 

Ip is the persistence length, and T is temperature (recall that Boltzmann's constant is unity, = 1); see sections 11 B 

for mode details. In particular, recall that for the Morse potential discussed around Eq. (26) the concrete formula for 

A reads A = l^-irp , where f ~ O.OleV, Ip ~ Inm and a Ip ~2 [7]. Taking room temperatures for T we get that A ~ 1. 
p 

Analogously, the dimensionless parameter /x enters the strand-surface attractive potential as ijV, such that the 

adsorbed phase of one single strand (that is without inter-strand interaction) exists for /i > 1; sec sections 11 A and 

lie for details. Note that /x has the same qualitative form (107), where u is the strength of the surface-strand 

potential. Recall that for the electrostatic surface-monomer attraction the concrete expression for is discussed in 

(22): n = fc^pj ^ r where is the Debye screening length (fc^^ = O.lnm at normal conditions), q is the monomer 

charge (aroimd one electron charge for a single-strand DNA), Ip is the persistence legth (around Inm for a single-strand 
DNA), and finally a is the charge density of the surface. Strongly charged surfaces have typically 1 — 10 electron 
charges per 1 nm^. At room temperatures /i ~ 1. 

In Fig. 1 the thcrmodynamical phases are confined by thick lines. ND, NA and DD refer, respectively, to the 
naturated-desorbed, naturated-adsorbed, and denaturated-desorbed phases. The meaning of these term should be 
self-explanatory, e.g., in the ND phase the two strands are localized next to each other, but they are far from the 
surface. 

First of all we see that there is no adsorbed and denaturated phase: as we have shown already in section IV, even 
small (but generic) inter-strand (inter-particle) attraction suffices to create a naturated state, provided that the two 
strands (particles) are adsorbed. Thus the rectangular region c in Fig. 1, which belongs to the naturated and adsorbed 
phase NA, refers to conditions where the overall binding is due to sufficiently strong attraction to the surface. 

The curved line going from (/i = 1, A = 0) to (/i = 0.25, A = 1) in Fig. 1 confines region a, where no overall binding 
(i.e., no-denaturation and no- adsorption) is possible according to the lower bounds obtained in the previous section. 

The region b, confined by two straight normal lines and the thick curve, refers to the the collective binding situation. 
It is seen that this region lies below both adsorption and denaturation thresholds. While we do not know the precise 
position of the thick curve confining the region b, we proved its existence in section V. 

Finally, the line separating NA (naturated-adsorbed) phase from ND (naturated-desorbed) phase extends mono- 
tonically to /z = 0.25 for A ^ oo. Please note that the monotonicity of this line is conjectured. Still this conjecture 
is, to our opinion, quite likely to be correct. 
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VIII. SUMMARY. 

The main purpose of this paper was in studying DNA denaturation in the presence of an adsorbing plane surface. 

As we argued in the introduction, there arc several relevant situations when the two processes, adsorption and 
denaturation, are encountered together. Taking into account the importance of these processes in the physics of DNA, 
as well as for DNA-based technologies, it is important to understand how specifically adsorption and denaturation 
interact with each other. 

Our two basic findings can be summarized as follows. First we saw that provided the two strands of DNA are (even 
weakly) adsorbed on the surface, there is no denaturation phase transition. There is only a smooth crossover from 
the naturated state to a (very) weakly bound state. Second we have shown that when the inter-strand attraction 
alone and the surface-strand attraction alone are too weak to create naturated and adsorbed state, respectively, their 
combined effect ("Borromean binding") can create such a naturated and adsorbed state. 

The results were displayed on a simple model of two coupled homopolymers (strands) interacting with the plane 
surface. The volume interaction within each homopolymer can be accounted for via renormalizing the persistence 
length; see section II F. Many realistic features of DNA are thereby put aside; see the beginning of section II. We 
plan to investigate some of them elsewhere. Another interesting subject is to study the DNA adsorption on a curved 
surface [28, 29]. 

We, nevertheless, hope that the basic qualitative aspects of the presented problem are caught adequately, and that 
the presented results increase our understanding of DNA physics. 
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APPENDIX A: DERIVATION OF THE SCHRODINGER EQUATION FROM TRANSFER MATRIX. 

Imposing the periodic boundary conditions, the partition function (3) can be written as 

Z = Trr^ = J dndf2T^{ri,f2;ri,r2), (Al) 



(A2) 



where T is the transfer operator parametrized with two continuous indices: 
T(fi,f2;fi',f2) = exp 
Thus in the thermodynamic limit N ^ oo: 

Z = A^, (A3) 

where A is the largest eigenvalue of T. 

For simplicity reasons, the subsequent discussion will be done in terms of a transfer matrix, which depends on a 
two scalar variables z' and z. The extension to the more general case (A2) is straightforward. 
Write the eigenvalue equation for the right eigenvector as 

J dz' e-'5^(^)-^(^-^')V(^') = e-^f'i^{z), (A4) 

where e"''-^ and ip are, respectively, eigenvalue and eigenvector. It is seen from (A3) that Nf is the free energy of the 
model in the thermodynamic limit N ^ 1 provided that there is a gap between the largest eigenvalue A and the one 

but largest eigenvalue. 
One now assumes that 

(3KD^ > 1, (A5) 

where D is the characteristic length of ^{z). Since ip^{z) is the density of monomers, we see that D quantifies the 
thickness of the adsorbed layer. Recalling (4) we can write condition (A5) as 

D^lp, (A6) 

i.e., the thickness is much larger than the persistence length. 

Under condition (A5) the dominant part of the integration in (A4) is z f» z' . With this in mind we expand V(-2') 
in (A4) as 

i,{z') = ^(z) + (^ - z')^\z) + ^i^£^V"(^) + (A7) 
and substitute this expansion into ( A4) . The outcome is 



2^ ^-m^) ( 1 + 1 - ) v,(^) = e-^^^{z). (A8) 



The corrections to this equation are of order Q( j^5^2£)4 ) = 0(;^). 
Eq. (A8) can be re-written as 

For weakly-bound states 

\V{z) - f\ « 1, (AlO) 
for those z, where |V'(-2)| is sufficiently far from zero. Thus in (A9) we can expand 

g/3(V(^)-7)_;^^^(y(^)_ J) (All) 

and get the Schrodinger equation: 

l^+f3^KV{z)^ ^{z) = (^/32k/+ ^ln|^) ^{z) ^ E^l^{z). (A12) 

The ground-state energy E of this Schrodinger equation relates to the free energy / of the original polymer problem. 
For weakly-bound states E is negative and close to zero. 
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APPENDIX B: SOLUTION OF SCHRODINGER EQUATION WITH THE DELTA-SHELL POTENTIAL. 



1. Discrete spectrum. 



Here we outline bound-state solutions of a one- dimensional Schrodinger equation 



-^'^"{^) + {y{^)-Emx) = Q (Bi) 



with the attractive delta-shell potential 



and with the infinite wall at a; = 0: 

?A(0) = 0. (B3) 

Here A > is the dimcnsionlcss strength of the potential, while xq is the radius of attraction, m is the particle mass. 
In contrast to the main text, hero wo do not put m = 1. 

Due to boundary condition (B3) the considered problem is equivalent to the corresponding three-dimensional 
Schrodinger problem with centrally-symmetric potential. 

Let us rewrite (Bl) as 



i)"{x)-k'^i){x) = -^di^x-xoji}, k = y/2m\E\ > 0. (B4) 

Xq 

For X ^ Xo (B4) is a free wave-equation. Its solution for x < xq and x> xq are found from the boundary conditions 
= 0) = and il}{x — > oo) = 0, respectively. Thus the overall solution is obtained as 

V'(a;) = M-""'^ sinh(fca;<) e-'=^> , (B5) 

a;< = min (a;, a:o)- x> = max(a;, xq), (B6) 
where N is the normalization constant determined via dxtp^{x) = 1, 



A/- e r ■ w^, X ^, 1 sinh (kxo) _2t^„ 

= — -- — smh(2A:a;o) - 2kxo] H — — e ^^'^^ . 

4k 2k 



(B7) 



Substituting (B5) into (B4) we get an equation for the energy of the single bound-state: 

kxQ = II sinh(fca;o) e''"'"'. (B8) 
The critical strength of the potential is seen to be 

11 = 1, (B9) 

because for small kxo Eq. (B8) gives 

kxo = (BIO) 
M 

Note the following form of ^^{x) for small values of kx^: 



Xq 

For large values of A the bound state energy increases as 

2kxo = II. (B12) 



^ We should like to clarify the physical meaning of studying the delta-shell potential (B2). First of all it should be clear that the weak- 
potential condition (All) docs not (formally) hold for the strongly singular potential (B2). Thus the transition from the transfer-matrix 
equation to the Schrodinger equation is foTrnally not legitimate. Nevertheless, there is a clear reason for studying the potential (B2) in 
the context of polymer physics, since it is known that the physics of weakly-bound quantum particles in a short-range binding potential 
does not depend on details of this potential [23] . So once the conditions for going from the transfer-matrix equation to the Schrodinger 
equation are satisfied for some short-range potential, one can employ the singular potential (B2) for modeling features of weaJfly-bound 
particles in that potential. This is in fact the standard idea of using singular potentials. 
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2. Continuous spectrum. 

For studying the continuous (positive-energy) spectrum of Eq. (Bl), we re- write it as 



^"(x) +n'^ij(x) = -—d(x- xo)tp, n = V2mE > 0. (B13) 



The solution is found as 



V2 

^{nx, n) = —= r- sin(na;<) sin(na;> + (5(n) ) (B14) 

VTT sin(na;o) 

where 5(n) is the phase-shift to be determined below, and where M is the normalization constant determined via 
orthogonalization on the n-scale: 



/ da;^/^(nx,n) ■i/'(na;,n') = (5(n — n'). 



(B15) 



This normaHzation can be checked via the large-a; behavior of 'ip{x,n) [23]. Note that for ^{nx,n) there are two 
types of dependence on the wave-number n: as a prefactor for the argument and as a parameter entering into the 
normalization and the phase-shift. 

For the phase-shift S{n) we get: 



that for n — > reduces to 



Thus 



and for tpinx, 0) we have 



n xq sin 5(n) = sin(n a;o) sin(n xq + 5{n) ), (B16) 



l_l = !^cot(<5(n)). (B17) 



5(0) = 0, (B18) 



V'(na:,0) = ^^"^^ sin(na;>). (B19) 

\/-KXq 



APPENDIX C 

While the finitcness of the integral (71) is obvious (because the integral dn/(l-|-n^) is already convergent), the 
convergence of the integral in (72) is less trivial. Estimating from (B14, B19) 

[2 

(j){nx,Q) = \ l — sin(na;), (CI) 

we get for the integral in (72): 

4 dnidr7.2 nlnl _ 1 /"^'^ dandn sin^ 2a 



TT^Jo \ + nl+nl [{nl-nlY + 2{nl + nl) + lf ^t^^ Jq Jo \ + n^ ^ cos2 2a + 2n2 + 1]' 

r'^ (Jq, (;Jj^ ^2 ^,jj^2 ^ 



47r2 7o 7o 5 + "- [n^cos^ Q + 2n+ 1]^ 



(C2) 
(C3) 



The integral over n in (C3) is convergent and produces an integrable logarithmic singularity ~ In cos^ a. Thus the 
double integral in (C3) is finite. 



